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1. Introduction and notations 



1.1. Aim and plan of the paper 

Techniques of integration over large matrices are important in several contexts of physics 
- from QCD M to quantum gravity [0, from disordered systems to mesoscopic physics 
|J- and of mathematics - enumerative combinatorics J|], integrable systems ]5],|6|,[7]] , free 
probability theory 0, statistics 0, etc. In most applications, a central role is played by 
the statistics of the eigenvalues of the random matrices - spectrum and correlations - or 
more generally of their invariants TrA p . This is justified once "angular" variables have 
been integrated over, like those appearing in TrAB. 
It is in this general context that the integral 



I=[ DU expTr AUBU* 

Ju(N) 



A and B hermitian, was studied more than twenty years ago, and an exact expression was 
derived Soon after, it was realised that this result had been obtained long before by 



Harish-Chandra |TT| as a corollary of a more general problem. The purpose of this paper is 
to return to this integral (sometimes called the Harish-Chandra-Itzykson-Zuber integral 
in the physics literature), to review the known facts and to present some of its known 
extensions, before turning to what remains a challenge: to find a good and systematic 
description of the large iV limit of its logarithm. We shall report on some recent progress 
made on the latter issue. 

Our paper is organised as follows. In the rest of sec. 1, we introduce notations and 
some basic results. The derivation of the expression of / is then reviewed in sec. 2, using 
various methods: heat equation, character expansion and Duistermaat-Heckman theorem. 
Sec. 3 discusses briefly extensions in various directions, in particular the case of rectangular 
matrices. Connection with integrable hierarchies is then presented in sec. 4, with special 
attention devoted to its dispersionless limit and what can be learnt from it. The resulting 
expressions for the large N limit of log / are then confronted in sec. 5 to those obtained 
by (what we believe to be) a novel diagrammatic expansion of log/ and by a purely 
combinatorial analysis [|12] . Finally sec. 6 contains a summary of results and tables. 
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1.2. Notations 

Let A and B be two N x N matrices. We shall assume that they are Hermitean, even 
though many properties that we shall derive do not require it. The subject of study is the 
integral 

I(A, B; s) = / DU exp — TrAUBU j (1.1) 

Ju(N) s 

where DU is the Haar measure on the unitary group U(N), normalised to J DU = 1; 
and the large N limit (in a sense defined below) of its logarithm ™ log I (A, B; s). At the 
possible price of a redefinition of U, one may always assume that A and B are diagonal, 
A = diag (ai)i<i<N, B = diag (i>i)i<i<jv- Clearly, the real parameter s could be scaled 
away. We find convenient to keep it as an indicator of the homogeneity in the a's and 6's. 

To any order of their 1/s expansions, both I(A, B\ s) and its logarithm are completely 
symmetric polynomials in the eigenvalues a^ and in the bi independently, as we shall see in 
the next section. We want to express them in terms of the elementary symmetric functions 
of the a's and of the o's: 

N N 
i=l i=l 

For finite N, only a finite number of these functions are independent, but this constraint 
disappears as N — > oo . By large N limit we therefore mean that we consider sequences 
of matrices A and B of size iV such that the moments of their spectral distributions 9 P 
and P converge as iV —>■ oo (see [|13| for some rigorous results on this limit). By abuse of 
notation we still denote by A and B such objets, and write 

F{A,B;s)= lim -L log I{A, B; s) . (1.3) 

We shall use various combinations of the symmetric functions 9 P and 9 P . In particular, 
for a h n, i.e. a a partition of n = J2 p P a pi which we also write a = [l ai ■ ■ -n a ™], we 
define 

tr Q A:=(-TrAj • • ■ (-Tr^j = J[ 6*. (1.4) 

p=i 

We shall also make use of the characters of the irreducible holomorphic representations 
of the linear group GL(iV), labelled by Young diagrams A with rows of lengths < Ai < 
A2 < • • • < Aat. They read 

XA = AW (L5) 
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in terms of Vandermonde determinants 

A(a)= J] (a l -a,) = det(ar 1 ) (1.6) 

l<j<i<N 

and their generalizations 

A A (o) = det(o^ +i_1 ) . (1.7) 
Frobenius formula |T^] relates these sets of symmetric functions: if a h n 

N^tr a A= J2 Xx(A)x\(a) (1.8) 

A 

| A | =71 

where the sum runs over all Young tableaux A with |A| = n boxes, and X\(a) denotes the 
character of the symmetric group & n , for the representation labelled by A and for the class 
labelled by a. 

The expression of integrals and differential operators over hermitian matrices in terms 
of their eigenvalues involves a Jacobian, and problems of normalisation appear. We discuss 
these questions shortly. 

If M = UAW, A diagonal, U unitary, we have dM = UdAW + [dX,M], where 
dX := dUU ] is antihermitian. Then Tr(dM) 2 = Ejdaf + 2 J2i<j \dXij\ 2 \a>i - aj\ 2 
defines the metric tensor g a p in the coordinates £ Q = (a^X^j). This determines 
first the measure DM = y/dBtgU^ a = 2 N{ ^ N - 1 ^ 2 ]\ i dM u ]\ i<3 d^eM i3 ;d3mMy = 
2 iV(7V-i)/2 A 2( a ) Yl den n dX tj = CA 2 (a) n< dcii DU. The constant C is fixed by comput- 
ing in two different ways the integral of a U(N) invariant function of M, for example a 
Gaussian 



e-^ TrM2 C f^r da, 



1 = / DM- = - / TT -^e-i< TT( 



di - a 3j 



The latter integral equals Ylp=i thus 

C= [ l N , , DM= ^ A 2 (a)nda,D^. (1.9) 

n p =ip! n p =ip! i= i 

From the metric above, one also computes the Laplacian 

1 d _ 9 1 9 T-r, 4 , 9 



= AM" 1 £(£)*(«) + * 



(1.10) 



where the last equality results from the vanishing of ^2 k [d ak , [9 afc) A(a)]] = 0, as this 
completely antisymmetric function of the a's is a polynomial of degree iV(iV — l)/2 — 2. 
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2. The exact expression of I(A, B; s) 

Assume that all the eigenvalues of A and B are distinct. One finds [ |11~| , |1C)|| 



N-l 



detfeT^), 

-l)/2 V '1; 



/(A, B- s) = ( J] P } -) (N/s)-^-^ A{aW lf' J ~ N (2-1) 

Note that both the numerator and the denominator of the r.h.s. are completely antisym- 
metric functions of the a's and of the b's independently, and that the limit where some 
eigenvalues coalesce is well defined. 

This expression ( |2.1| ) may be obtained by several different routes. 

2.1. Heat equation 

For two hermitian N x N matrices M A and Mb, let us consider 

/ M \ N2/2 N 
K(M Al M B ;s) = — exp--Tr(M A -M B ) 2 . (2.2) 

\Ztts I is 

K{M Al Mb; s) satisfies the heat equation 

d 1 



(h 2 A Ma \K(M A ,M B ;s)=0 , (2.3) 

where Am is the Laplacian over M, together with the boundary condition that for t — > 0, 
K(M A , M B ;s)-+ 5(M a - M B ). 

The heat kernel K(M A , Mb', s) is invariant under the simultaneous adjoint action 
on M A and Mb by the same unitary matrix U. If we diagonalise M A = U A AU A and 
M B = U B BU B , K(M A: M B ; s) = K(A, UBW; s) where U = U A U B - Upon integration 

K(A,B;s):= j DU K(M Al UM B U ] ; s) = J DU K(A, UBU ] ;s) 



is again a solution of the heat equation ( |2.3| ) , but depends only on the eigenvalues of A 
(and hi of B). Using the explicit form of the Laplacian in terms of the eigenvalues a^, (see 
( OP ), K satisfies 



The product A(a)A(b)K(A, B; s) is an antisymmetric function of the a's and of the 6's, is a 
solution of the heat equation with the flat Laplacian, and satisfies the boundary conditions 
that for s — > 0, CA(a)A(b)K(A, B; s) — > jjj J2pee N ep Hi $( a i — bpi), with C the constant 
computed in ( |1.9|) . 

In physical terms it is the Green function of N independent free fermions, and it is 
thus given by the Slater determinant 



1 / N \ N/2 
CA(a)A(b)K(A, B;s) = — - det 



exp-^(a; - bj) 2 



(2.6) 



which is consistent with (12. If). 



2.2. Character expansion 

We may expand the exponential in (|1.1|) to get 

I(A, B;s) = jr (iV/ ^ n I DU(TtAUBU^) n . (2.7) 
n=o n ' J 

Frobenius formula ( |1.8| ) evaluated for the partition [l n ] gives for any matrix X of GL(iV) 

Tr n X = dxXx(X) , (2.8) 



A 
1*1 = 



where d\ = X\([l n ]) is the dimension of the A-representation of & n . Integration over the 
unitary group then yields 

f BU Xx (AUBUi) = XA^)XA(fl) = Xx(A)xx(B) 
J ' Xx(I) d x 

and a well known formula jI3| gives 



Using ( |1.5| ) and putting everything together we find 



N-l oo 

A(a) A(6) I(A, B; S ) = ( J] pi) £ n fA + » - IV Aa(q)Aa(6) (2 " 11( 

p=l n =0 * LLpK P~1~P )■ 



V- 

|A|=n 



Ar_1 W /"AT"/ ^ -0+1 

(n?o e n ( det^o^) (2.1m 



/,-: n<\/ ; ...1 -v </-] ^ 9 
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Eq. ( [2.11a ) is interesting in its own sake, while an extension of Binet-Cauchy theoremS 



enables one to resum ( |2.116p into 



N-l 



A(o) A(6) I{A, B; s) = ( J] pi) (N/s)- N(N - 1)/2 det (e^) (2.12) 



P =i 

which is precisely 

£.5. Duistermaat-Heckman theorem 



Let us compute the integral by the stationary phase method. The stationary points 
U of the "action" Tr AUBU^ satisfy Tr 6UU& [U BU£, A] = for arbitrary antihermitian 
5UUq, hence [UqBUq, A] = . For diagonal matrices A and B with distinct eigenvalues, 
this implies that UqBUq is diagonal and therefore that the saddle point Uq are permutation 
matrices. 

Gaussian fluctuations around the stationary point Uq = P may be computed by 
writing U = e x P, X antihermitian, and by integrating over X after expanding the action to 
second order. Summing over all stationary points thus gives the "one-loop approximation" 



to the integral 



I(A,B;s) = C £ e-H>^ /n^ e "^ El<j|X " |2(ai " aj)(bPl " bPj) 

D r~ (Cr v • ^ ■ 



Pee N " i<j 

■ 7rs ,N(N-l)/2 



p7^ n HiKjiai-ajXPpi-b 



t<j yui - Ujjyupi - u Pj/ 



N) A(o)A(6) 

- c ' (S) ^Wr (2 - 13) 

1 Recall P that if f(x) = J2 e > Q f&'i 

//i--- /e w det det b\ 3 = ^ ^ ^2 ep £p ' ^ N II a P* b P^ 

' ii>OP,P' i 

Jf\ £P F " II S ^ ^P-P'i^ 

' P,P' i li 

det f{aibj) 



o<e 1 <e 2 <--<tN ii>op,p> 
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which reproduces the previous result up to a constant C . The latter may be determined, 
for example by considering the s — > limit, and the result reproduces ( |1 . 1| ) . Thus the 
stationary phase approximation of the original integral ( |1 . lj ) or ( |2.4j) turns out to give the 
exact result! This well known empirical fact turned out to be a particular case of a general 
situation analysed later by Duistermaat and Heckman [ [L6[ : if a classical system has only 
periodic trajectories with the same period, the stationary phase (or saddle point) method 
is exact. 

In more precise mathematical terms, let M. be a 2n-dimensional symplectic manifold 
with symplectic form u, and suppose that it is invariant under a U(l) action. Let H be the 
Hamiltonian corresponding to this action (i.e. dH = iyu, V the vector field of infinitesimal 
£7(1) action). Assume also that the fixed point set (the critical points) is discrete. Then 
the theorem of Duistermaat-Heckman asserts that the stationary phase method is exact, 
i.e. that 



e itH 



^L] n e *i si gn (Hess(p c )) el tH(p c ) Vdetu(P c ) 
t J Y c v/|detHess(P c )| 



where the sum is over (isolated) critical points P c ; the phase involves the signature 
sign(Hess(P c )), i.e. the number of positive minus the number of negative eigenvalues, 
of the Hessian matrix Hessij = d 2 H/d^d^ 

The integral (|1.1|) satisfies the conditions of the above theorem. The integration runs 
over the orbit O = {M = UBU^} of B under the coadjoint action of U: this orbit, home- 
omorphic to the manifold U(N)/U(1) N , has the even dimension N(N — 1) and is in fact 
a symplectic manifold. On two tangent vectors Vi = [Xi, M], i = 1, 2, (JQ antihermitian) , 
tangent to O on M, the symplectic form reads u>(Vi, V2) = TrM[Xi,X2]. The Hamilto- 
nian H = TrAM defines a periodic flow M(t) = e lAt M (0)e~ lAt if all eigenvalues of A are 
relatively rational. As the latter configurations form a dense set among diagonal matrices 
A, this constraint may in fact be removed and this justifies a posteriori the stationary 
phase calculation above. For further details on the Duistermaat-Heckman theorem in the 



present context, the reader may consult also []TT| , P~8 



3. Generalizations and extensions 

3.1. Other groups 



As already mentioned, the integral ( |2.1|) appeared first in the work of Harish-Chandra JTT 



as a simple application to compact groups of a general discussion of invariant differential 
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operators on Lie algebras. Following Harish-Chandra, let G be a compact Lie group. 
Denote by Ad the adjoint action of G on itself and on its Lie algebra g, by (.,.) the 
invariant bilinear form on g, and by f) C g the Cartan subalgebra. If hi and h 2 £ f) 

A(hi)A(h 2 ) / dg exp(Ad(fir)/ii, h 2 ) = const e(ty) exp(«;(/ii), /12) , (3-1) 

where w is summed over the Weyl group W, e(w) = (—l) i ( w \ £(w) is the number of 
reflections generating u>, and A(h) = ria>o( a ' h) ' a P r °duct over the positive roots of g. 
In the case of U(N), if we take hi = iA, h 2 = iB, (|3.1| ) reduces to ( |2.1| ). 

This extension to general compact G may also be derived by constructing the heat kernel 
K(gi, Ad(g)g2] s) on G in terms of characters, by using Weyl's formulae for characters and 
by averaging K over G Jl9]. Let's sketch the derivation. One may always assume that 
9j = e i & i), j = and one writes 

K{g l: g 2 -s) =Y J dxXx{9i92 1 )^ sCx (3-2) 

A 

where C\ is the value of quadratic Casimir for the representation of weight A, d\ the 
dimension of the latter. Integration over the adjoint action gives 



dgK(g u Ad(g)g 2 ;s) = ]T XxMxxig^e-^ • (3.3) 

A 

Weyl's formula for the characters reads 

X\(e ) — n 7 n , (3.4j 

same tor A = U 

where p is the Weyl vector, sum of all fundamental weights. As C\ = |A + p| 2 — \p\ 2 , the 
r.h.s of ( p.3| ) is the exponential of a quadratic form in A. The summation over A may be 
extended from the Weyl chamber to the full weight lattice, Poisson summation formula is 
then used and after taking a limit of infinitesimal s, hi and h 2 , one is led to (|3.1|) . 

For completeness we also mention the generalisation of ( |1 . 1|) involving Grassmannian 
coordinates and the integration over the supergroup U(Ni\N 2 ), see [|0],^l[. The case of 



integration over a pseudounitary group U(Ni,N 2 ) has also been discussed, see [[22) and 
further references therein. 
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3.2. Rectangular matrices 
Consider the integral 

I (2) (A 1 B;s)= [ DU I DFexp — Tr(AUBV* +h.c.) (3.5) 

JU(N 2 ) Ju(N!) s 

where A is a complex JVi x N2 matrix, B a complex N2 x N\ matrix, and N = min(iVi, A2). 
Without loss of generality, one may assume that N\ > A2 and that the N2 x N2 matrices 
A^A and BB' are diagonal, with real non negative eigenvalues a^, bi respectively. We 
assume again that the a's on the one hand, the b's on the other, are all distinct, so that 
neither of the Vandermonde N2 x A^ determinants A(a), A(6), vanishes. Then the methods 
of heat equation or of character expansion presented in sec. 2 yield the following expression 

A(a)A(b)I^(A,B;s) = ^^1]^ , 

(3.6) 

with /^(z) the Bessel function /^(z) = Xl^Lo n!(n+t/)! (f ) U+U ■ This expression has been 
obtained for N\ = N2 in [pi, 15] and in the general case in [23|| . 

Integrals (1.1) and (3.5) are the cases i^T = 1,2 of an infinite set of unitary integrals 
which are exactly calculable: 

K A" 

/W(A fc ,S fc ; S )= / \\DUkexp-Y' Tr A k U k+1 B k Ul 



1 fc=l 



'u k eu(N k ) fe= 

= ^^^y (3.7) 

A(o)A(6) 

(index fc is cyclic modulo X) where A 7 = minfc(A^fc), and are A 7 ^ x N k+ i ma- 
trices, 1 < k < K; the generalized hypergeometric series is given by (j)(x) = 
^2^Lq xU / Y[ k =i( n + -^fc — A 7 ")!; and, assuming A 7 ! = A 7 , the cij (resp. bi) are the A 7 (dis- 
tinct) eigenvalues of AiA 2 . . . Ajc (resp. Bk ■ ■ ■ B 2 Bi). 

In what follows, we shall concentrate on integrals ( |1 . 1| ) and (|3.5|) ; however the analysis 
applies equally well to the more general integral ( |3.7|) , see in particular for a discussion 
of hypergeometric tau-functions of Toda lattice. 

3.3. Correlation functions 

Returning to the unitary group and the integral (|1.1|), it is also natural to consider the 

correlation functions associated with it, i.e. to compute the integrals 
f N 

/ DU U ilh ■ ■ ■ U lmJm Ut i£i ■ ■ ■ Xj\ mim exp -Tr AUBrf . (3.8) 

JU(N) s 

Partial results have been obtained in [24| and in ||25|| . We still lack explicit and general 



expressions for these correlation functions. 
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4. Connection with integrable hierarchies 



The integral ( |1 . 1|) turns out to provide a non-trivial solution of the two-dimensional Toda 
lattice hierarchy J26| . This stems from the following observation: define 



det (eH°* 6 i) . </v 
T7V " A(a) A(6) 

where h = s/N. Comparing with Eq. ( |2.1|) , we see that r^r = h~ N ( N ~ 1 ^ 2 n^Lo 1 ^') -1 ^- 
Then the following formula holds: (see also |26|,[Z7| for a two-matrix-model formulation) 



TN = d et(//^^ M VeU^ 1 ^ + E 55I V' + «- 1 «- 1 ) N ) (42) 
\J J 2niu 2niv j 

V J / 0<i,j<N-l 

where the integration contours are small enough circles around the origin, and with the 
traditional notations: 

N N 

t q = K-J2a$ = -0 q t q = h-Tb^ = -e q q>l (4.3) 
q q q q 

i=l ' i=l 

Formula ( fl.2|) can be easily proved by noting that e h ^- J i tqU = n^Li(l ~~ uai)~ l , 
e i E, = n « i(1 _ „,,)-!, and expanding the contours t0 catch the potes at „ = a~\ 

v = b~ x . This makes tn a tau function of the 2D Toda lattice hierarchy, as we shall discuss 



% 

now. 



4-1. Biorthogonal polynomials and 2D Toda lattice hierarchy 

Noting that the parameter h = s/N can always be scaled away we set h = 1 throughout 
this section. We take Eq. (|4.2|) to be the definition of tn as a function of the two infinite 
sets of times (t q , i q ), q> 1. We also set To = 1. 

Formula ( f4.2| ) suggests to introduce a non-degenerate bilinear form on the space of 
polynomials by 

(q\p) = <p<p p(u)q(v)e^i>^ q q (4.4) 

x 1 ' J J 2mu 2-kiv 

and normalized biorthogonal polynomials q n (v) and p n (u) with respect to the bilinear form 
above, that is polynomials of the form p n (u) = h~ x u n H — • and q n (v) = h~ x v n + • ■ ■, such 
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that (q m YPn) = fimn for all to, n > 1. One can now replace monomials in Eq. (4J2) with 
biorthogonal polynomials and obtain immediately 

N-l 

tn = n h i ■ ( 4 - 5 ) 

i=0 

Next we introduce the matrices of multiplication by u and v in the basis of biorthogonal 
polynomials: 

Lmn = (q m \ U \p n ) L mn = (q m I V \p n ) TO, 71 > 1 . (4.6) 

Note that 

Lmn = TO > 71+1, L mn = 7i>TO+l. (4.7) 



A standard calculation []26| leads to the following evolution equations for L and L 
with respect to variations of the t q , t q : 

g = g = [(*)-.*] (««) 

S^"^ 1 ] 9*^ = [(i')-.i] ■ (4.8*) 

Here (-)± denotes the lower / upper triangular part plus one half of the diagonal part. 
Eqs. ( |4.8|) are the standard form of the two-dimensional Toda lattice hierarchy |28| (up to 
a choice of sign of the t q ). 

It is often more convenient to write these equations as quadratic equations in the set 
of (tjv)- These are the "bilinear" Hirota equations. They are obtained by picking two sets 
of times (x q ,x q ) and (y q ,y q ) and writing in two different ways 

du dv , . . . V s-w'+y y q v q +u~ 1 v- 1 

2mu zttiv 

where p n ,x,x is the right biorthogonal polynomial associated to times (x q ,x q ), whereas 
Qm, y ,y is the left biorthogonal polynomial associated to times (y q ,y q ). One finds: 



p-u n - m r n (x q - -u-\x q )r m+l {y q + l u -^y q )^>^ ^ 
^v m - n r m (y q ,y q - -v^)r n+1 (x q ,x q + ^)e^>^"-^ )v9 . (4.9) 
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Equivalently, equation (O) can be derived directly from the original expression ( |4.1|) using 



the set of determinant identities 

n+l m+1 

det (^ -a iy B) det(A' + a u B') = ^ det(A', B' - b' 3 ) det(A, B + 6j) (4.10) 

i=i i=i 

where A = {ai, . . . , a n+ i}, S = {&i, . . . , 6 n }, A' = {a' l5 . . . , a' m }, 5' = {b' x , . . . , b' m+1 }, and 
a symbolic notation det(-, •) is used for determinants of the form det(exp( j-a;6j)). 

Expanding Eq. ( ^.9| ) in powers of x q — y q and x q — y q , expanding in power series in u, 
v and performing the integration results in an infinite set of partial differential equations 
satisfied by the (tn). 

Example: Expand to first order in x\ — yi, and set m = n + 1. The result is: 

r n+1 T n - 1 = r n ddr n - dT n dr n Vn > 1 (4.11) 

(with d = d/dti, 8 = d/dtx) which is a form of the Toda lattice equation. It is of course 
also the Desnanot-Jacobi determinant identity applied to Eq. ( |4.2| ). 

Finally, the matrices L and L satisfy an additional relation, the so-called string equa- 



tion, which takes the form [26 



[L-\L~ 1 ] = 1. (4.12) 
4-2. Large N limit as dispersionless limit 

We now restore the parameter h, which is required for the large iV limit. Indeed, as 
iV — > oo, h must be sent to zero, keeping s = hN fixed. We define 

f(t q ,t q ;s)= lim h 2 \ogT N (t q J q ) (4.13) 

N— >oo 

where TAr(t g , t q ) is defined by Eq. (|4.2|) . In a region of the space of parameters (t q , t q ) which 
includes a neighborhood of the origin t q = t q = 0, f(t q , t q ; s) is governed by "dispersionless" 
equations we shall describe now. f(t q ,i q ; s) is called in this context the dispersionless tau- 
function. 

First, comparing with the expression (|1.3|) of the free energy F(9 q s~ q ^ 2 ,9 q s~ q ^ 2 ) 
(where we have explicitly stated the dependence of F on the 9 q and the 9 q and scaled 
away the parameter s), using Eq. ( |2.1| ) and definitions ([Of), we see that 



f(t q J q ;s) = ^s 2 \ogs + ^s 2 + s 2 F(qt q s- q / 2 -\qi q s- q / 2 - 1 ) . (4.14) 
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This is the scaling form of the dispersionless tau function.! 

The dispersionless Toda hierarchy is a classical limit of the Toda hierarchy: it is 
obtained by replacing commutators with Poisson brackets defined by 

{g(z,s),h(z,s)} = zf z f s -z^f s (4.15) 

Here z is the classical analogue of the shift operator = 6{j+i, which justifies that 
{logz, s} = 1. 

The Lax Eqs. ( |4.8| ) become 

j£ = -{(0+.4 W q =m ~' i} (4 - 16a) 

§- q = -m + A w q = {{iq) - J} {4 - m) 

where (-)± now refers to positive and negative parts of Laurent expansion in z, and £(z, s, t) 
and £(z, s,t) have the following z dependence: 



oc 

r z + ^\ k z~ k (4.17a) 

fc=0 

oo 

l=rz~ 1 + J2^kZ k (4.176) 

fc=0 

related to the structure ( |4.7| ) of L and L. 

These equations imply that the differential forms d<p(£, s, t), d(p(£, s, t) are closed: 

dp = md£/£ + logzds + J2( £q )+ dt q + X^ 9 )- d *« ( 4 - 18a ^ 
dip = rhd£/£ + logzds + ^(f)+di 9 + ^(^)_dt, (4.186) 

<?>1 ?>1 

where m and m are Orlov-Shulman functions, which can be characterized by 

m = J2qt q £ q + s + J2^ rQ ( 4 - 19a ) 



q>l q>l q 



dt q 
df 



m = J2qi q £ q + s + J2dT i ~ q ( 4 - 1% ) 

q>l q>l q 



2 In fact, one can get rid of one more parameter since A and B can be scaled independently, 
but we choose not to do so for symmetry reasons. 
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and satisfy the "dressed" Poisson bracket relations 

{log£,m} = l {log£,m} = l. (4.20) 

In the present case, we have the following constraints, which determine uniquely the 
solution of the dispersionless Toda hierarchy: 

m = m = (ii)- 1 . (4.21) 

Eqs. ( f4.21|) are directly related (via Eqs. (|4.20| )) to the classical limit of the string equa- 
tion (ETLJ), i.e. {i- 1 ,!' 1 } = 1. 

Let us call a = £~ x and b = The fact that m = fa implies^ that (p(a) and (p(b) are 
related by Legendre transform, or that their derivatives b(a) = -^(p(a) and a(b) = -^<p(b) 
are functional inverses of each other (the latter fact was derived by ad hoc methods in 
and 0). 



4-3. Application of the dispersionless formalism 

In what follows we set s = 1, so that 9 q = qt q . In order to explore the structure of the 
function f(t q ,i q ;s), we now assume that only a finite number of t q and i q is non-zero. 
Note that this cannot happen if the eigenvalues are real; however, here we are interested 
in properties of the integral ( |1 . 1| ) as a formal power series in the t q and therefore we do 
not worry about the actual support of the eigenvalues. Then, according to Eqs. (|4.17|) and 
( |4.19|) , the Laurent expansion of a{z) and b(z) is finite, of the form: 



n+l n+1 

a = J2 a Q z ~ q b=^f3 q z q (4.22) 

5 = 1 q=l 

where n = max{q\t q ^ 0}, n = max{q\i q ^ 0}, ol\ = (3\ = 1/r, a s+ i = nt^r n+1 , 
(3 n+ i = nt n r n+1 (r is defined by Eqs. (|4.17|) ). Then it is easy to show that a and b satisfy 
an algebraic equation, some coefficients of which can be worked out explicitly: 

b n+l a n+l _ b n a n _ b n( a n-l ei + , . . + £j _ fl n + . . . + £_) +lower Qrder termg = Q _ ( 42 3) 



Plugging Eq. ( |4.22| ) into Eq. (|4.23| ) leads to algebraic equations for the coefficients a q , (3 q 



as a function of the 9 q , 9 q . 

3 The constraint m = m is true for a more general class of solutions of Toda, see for example 



sec. 4.4 and J29]. 
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Examples: 

• If O q = 5 q x9i, one can set 6\ = 1 by homogeneity. In this case Eq. ( |4.23|) is quadratic 
in b. In terms of b and I = 1/a it reads 



b z - b(£ + d x t + ■■■ + 6 n £ n+L ) - £P(£) = (4.24) 

where P is a polynomial with -P(O) = 1. One way to determine P is to note that Eq. ( |4.22| ) 
provides a rational parameterization of a and 6, so that the resulting curve (a, b) has genus 
zero, which forces the discriminant of Eq. (|4.24|) to be of the form £(1 + £/ (4i(j 2 ))Q(£) 2 
where Q is a polynomial of degree n with Q(0) = 2, and ifi is a constant. This yields the 
expression 

&(*) = \(* + eif + ■ ■ • + ^ r+1 + 4^) QW ) ' (4 ' 25) 

Considering Eq. (|4.19a|) as an asymptotic expansion of m = b/£ as £ — > oo fixes Q: Q is 
the polynomial part of 2^(1 + 6>i£ H h 6> n £ n )/A/l + 4^ 2 /^- Imposing Q(0) = 2 leads to 

^ = 1 + E(-1) 9+1 ^^ 29+1 (4.26) 

Note that ip = r 2 = d 2 F/d6 1 d6 1 = exp(d 2 F/ds 2 ) (the latter equality being the disper- 
sionless limit of Toda Eq. ( fl.ll| )). 

By Lagrange inversion, Eq. ( [4 . 2 6|) allows to get the exact expansion of if) as well as of 
F in powers of the 6 q . One finds, restoring the 6\ dependence: 

• If 6> g = 5 qn 6 n , 9 q = 5 qn 9 n (with the same n), one can show using the formalism 
above that ip = d 2 F jdQ\dQ\ satisfies the equation: 

^=l + (n+l)Mnf +2 (4.28) 

so that 

g ((n + l)fc + l)!fc! g ((n + l)fc)!fc! 

(4.29) 

For example, for n = 2, F = ±(6> 2 # 2 ) + f (6> 2 # 2 ) 2 + §(6> 2 # 2 ) 3 + • - • 
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4-4- Case of rectangular matrices 

All of the formalism above applies equally well to the integral ( |3.5| ) over rectangular ma- 
trices. The latter has a "diagonal" character expansion i.e. 

r (2) - const I™(A B- s) - V^"I A I - Aa(g) A ^ 6 ) 

r N -const { A,B,s) Up (X p + p - l)l(X p + p - 1 + u)l A(a) A(6) 

(4.30) 

where H = s/N, N = min(iV"i, N 2 ), v = \N X - N 2 \. 

The form of Eq. Q4.30| ) alone proves that the ("rjy ) form a tau function of 2D Toda 
lattice (for fixed v). There is also a determinant formula similar to Eq. (|4.2|) which we 
shall not write. 

The large N limit is taken keeping hN = s and fiv = £ fixed. The times are 
t q = h-^Tr(AA^) q , t q = h^Tt(BB^) q . Biorthogonal polynomials and the dispersionless 
equations can be obtained similarly as above. 

We have thus constructed a whole family g^(t qi i q ; s) = limjv^oo fo 2 logrj^ of disper- 
sionless tau functions. The diagonal character expansion ( f4.30| ) implies that they belong 
to the class of solutions that satisfy m = rh. After some calculations, one finds that the 
constraint analogous to Eq. (|4.21|) is 

m(m + = m(m + = (U)' 1 (4.31) 

Note that if £ = i.e. N\ = N 2 , the constraint (|4.31| ) becomes m = rh = (££) -1 / 2 , 
which suggests to redefine I = 1% , I = l 2 ■ The new functions £ 2 , l 2 satisfy the string 
equation {£ 2 , £ 2 } = 1 of the integral (|1.1|), and of course the same Toda equations ( |4.16| ) 



for even times, with the replacement t 2q — > t q , t 2q — > i q . In fact, one can check that we 
recover in this case the dispersionless tau function of the usual integral ( |1 . 1| ) with odd 
moments equal to zero: g (u q ,u q ;s) = 2f(t 2q = u q ,t 2q+1 = 0,i 2q = u q ,t 2q+ i = 0;s). 

5. Diagrammatic expansion and combinatorics 

5.1. A Feynman diagram expansion of F 

To develop a Feynman diagram expansion of (|1 . 1] ) , we first trade the integration over the 
unitary group for an integration over N x N complex matrices X by writing 



^Tr AUBU 
e s 



f = (*) ~ n2 [ DXDX^e- NT * xxi exp ^=Tr(AUX^ + XBU^) (5.1) 
AT J y/s 
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Thus, 



where 



I(A,B;s) = const j DX DX^- nt ^ xx ^ +n2w ( axbx " » 



(5.2) 



^N 2 W(X 1 X 2 ;s) 



N , 
DUexp — TrOJXt + U^X 2 



This integral is known exactly in this large iV limit 
as mentionned in sec. 1, 



(5.3) 

: with the same abuse of notations 



W(Y;s) 



oo 

E 

n=l 



tr„F 



a hn 



{2n + J2a p -3)\ 
(2n)! 



n 



~(2p)! 
pl(p — 1)! 



(5.4a) 
(5.46) 



and W a is an integer, as follows from the recursion formulae discussed in |^ . 

The form fl5.2j) is adequate to develop a diagrammatic expansion a la Feynman of 
F(A,B;s) (see fl3~3| , |34|| for reviews). Recall that double lines are conveniently introduced 
to encode the conservation of indices []35| . The inverse of the first term in the exponential 
of (|5.2| ) yields the propagator, (X^X^) = \ \ [ = j^5u5jk, while each term in W, 
i.e. each monomial tr a (AXBX')/ Y[ p (&p\p otp ) gives rise to a multi-vertex of type a, 
(see Fig. 1), which comes with a weig ht iv 2 -E Q p W a times products of a's and 6's and 
Kronecker symbols expressing the conservation of indices. When these Kronecker symbols 
are "contracted" , they leave sums over closed circuits (or faces) of powers of a's or 6's. A 



face of side p thus contributes a factor N9 P , resp. N9 q with the notation of (|1.2| ). To put it 
another way, the Feynman diagrams may be regarded as bicolourable, with faces carrying 
alternatingly a or b "colour" . 

As is well known, only connected diagrams contribute to the free energy log /. In the 
discussion of this connectivity, each multi-vertex just introduced must be regarded as a 
connected object, and it is useful to keep track of this fact by drawing a tree of dotted 
lines which connect the various traces which compose it, see Fig. l.b. 



Ik 




Fig. 1: (a) a vertex of type [2]; (b) a 
multi- vertex of type [l 2 2 1 3 1 ]. 



(a) 



(b) 
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If a diagram contributing to log / has P propagators, V a vertices of type a and L 
loops of indices, and builds a surface of total genus g with c connected components, it 
carries a power of N equal to 

#( jV) = -P + V a (2 - Yl a p) + L 

a p 
ct p a p 

= 2 - 2g + 2 [c - 1 - J2 V*(%2 a P - X )) ^ 2 " 2 3 > 

a p 

where the last inequality expresses that to separate the diagram into c connected compo- 
nents, one must cut at most its Yl a ^q(S p a p ~~ 1) dotted lines. 

This simple counting (see the second ref ||32|| ) shows that the leading 0(N 2 ) terms 
contributing to F are obtained as the sum of planar (i.e. genus g = 0) and minimally 
connected Feynman diagrams, such that cutting any dotted line makes them disconnected. 
They are thus trees in these dotted lines. 

These Feynman rules, supplemented by the usual prescriptions for the symmetry factor 
of each diagram (equal to the inverse of the order of its automorphism group) , are what is 
needed to get a systematic expansion of F in powers of the 9 and 9. 

f = e ixjj^r n w « n h n v ^ 

minimally connected I I vertices a— faces b— faces 

diagrams Q 

Note that these rules imply that the sign attached to each diagram (coming from the 
signs of the W a ) is also (-1) L = (-l)P ower of d + P ower of ~ e . 

Examples: 

W(XX^ ; s) contains a term linear in XX^ (with coefficient 1/s), which may be in- 
serted in arbitrary number on any propagator without changing the topology of the dia- 
gram: each individual insertion (which raises the corresponding power of a and b by one 
unit) is depicted by a cross on the propagator. 

Terms in F linear in 9 and 9 come solely from these insertions, which yields 
SpLi p"^- \ jf ( a single loop with p insertions has a symmetry factor equal to p). 
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Terms in F of the form (Qi9\/ s) p come entirely from a multivertex of type iT^ V ]XX\ 



whence a contribution 2 P ^th^t W „ • More generally, the term in F of the 



form rip^p) ^! 1 comes solely from the multivertex of type a, a h n, in which the con- 
traction of b lines into "petals" is unique and determines that of a lines. The weight W a 
together with the symmetry factor n p P apQf p' reproduces the result of Ql27| ). 

Likewise, it is easy to find the diagrams contributing to the {O2O2/ s 2 ) 2 and (O2O2/S 2 ) 3 

terms in F ^ /^^^ which give (± + (±) 2 )(6> 2 # 2 /s 2 ) 2 = j(0 2 8 2 /s 2 ) 2 and 

■ 




■>"■> \ 9 f 02^2 



. # 

in agreement with the expansion of the end of sec. 4. Conversely, one sees how effective 
the methods of sec. 4 are to resum classes of Feynman diagrams. Indeed it is not obvious 
how to derive directly from the diagrammatic expansion that -^(#2^2)5 or more generally 
F{6 n 9 n ) and i/j(9 n 9 n ), have the simple form given by Eq. (|4.29|) . In particular it would be 
interesting to find a direct combinatorial proof of the equation ( |4.28|) satisfied by ip. Its 
form suggests a possible connection with decorated rooted trees, perhaps in the spirit of 
0. 

Remark: Actions of the form of Eq. ( |5.2j ) and their diagrammatic expansion are a 
generalization of the "dually weighted graphs" of ||37|| : the extra ingredient is that graphs 
are required to be bi colored and the black/ white faces are weighted separately. They 



are also a generalization of the bicolored diagrams of |38| - the faces of the latter being 
weighted only according to their color and not to their number of edges (this corresponds 
to the case where A and B are projectors). 

5.2. Combinatorics 



We now show that this expansion is in fact equivalent to the one recently obtained in |T^ . 



What is remarkable is that the methods are orthogonal: [T2\ is based on manipulations 
in the symmetric group, and in particular uses some results on the number of solutions of 
equations for permutations (which are however themselves related to planar constellations 
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||36|| ) . We now explain the results of [12 1 that are relevant here, and show their equivalence 



to our Feynman diagram expansion. 

First we explain how to associate to a (not necessarily connected) bicolored map g, a 
pair of permutations a(g) and r(g) of the set of edges. a(g) (resp. r(g)) is the permutation 
which to an edge associates the next edge obtained by clockwise rotation around the white 
(resp. black) vertex to which it is connected, see the example of Fig. 2a). This is a one-to- 
one correspondence. 

An important remark is that the permutation a(g)r(g) encodes the faces of the map. 
Indeed, cycles of a(g)r(g) are obtained by turning clockwise around each face and keeping 
only the edges going from a white to a black vertex. It is implied that each connected 
component has its own face at infinity (for which rotation must be made counterclockwise). 
Consequently, the lengths of the cycles of a(g)r(g) are precisely one half of the sizes of the 
faces. 



Fig. 2: a) A bicolored planar map and b) its dual, a = (1 2 11) (3 4 5) (6 7 10) (8 9), 
r = (1) (2 3 6) (4 5) (7 8) (9 11 10), err = (1 3 5 6 8 11) (2 10) (4) (7 9). 



Also, for p a permutation of the set of edges, define H p to be the partition of the set of 
edges into the orbits (cycles) of p. For two such partitions IT and IT', we say that II < IT iff 
every block of IT is included in a block of IT'; and define IT V IT' to be the smallest common 
majorant. Also, call #11 the number of blocks of EE. Finally, if n is the number of edges, 
define C p h n to be the partition of the integer n corresponding to i.e. the lengths of 
the cycles of p. 

Then the results of [[T2] (Theorems 2.12, 4.2) can be reformulated as follows: 




b) 



oc 



— n 






n=l 



#n CT +#n. r +#n (TT -n=2#(n (T vn r ) 
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The coefficient j(p, II) is only defined explicitly if II = n p , in which case ^(p, H p ) = Wc p 
where W a is given by Eq. ( |5.4fc|) . In general, we have the following expression for 7(p, IT): 

7 (p,n)= Yl II w ^ ( 5 - 7 ) 

n'>n p ,#(nvn')=i blocks i 
#n p -#n'=#n-i of n ' 

where the pi are the restrictions of p to each block i of the partition IT'. This expression 
looks somewhat complicated but it is very easy to interpret once we have related this 
formalism to the diagrammatic expansion of the previous section. 

In order to proceed with the equivalence, we first associate to the pair of permutations 
a and r a map according to the construction above. Since a and r are permutations of 
{l,...,n}, this produces a map with labelled edges from 1 to n. The quantity + 
$TI T + ftHfTT — n is simply the Euler-Poincare characteristic of the map and the condition 
on the summation simply imposes planarity of the map (or more precisely, of each of its 
connected components). 

Secondly, we unlabel the resulting map. That is, since the summand of Eq. ( |5.6| ) 
does not actually depend on the labelling, one can sum together maps which are only 
distinguished by the labelling of edges; by definition of the symmetry factor of a map 
(inverse of the order of the automorphism group, i.e. here number of permutations that 
commute with both a and r), the only modification this produces is to replace the 1/n! 
with the symmetry factor of the unlabelled map. Thus, we rewrite Eq. (|5.6|): 



s 

f ( A ' 5 ' s ) = E^^H^W,n ff ( 9 )Vn T(9) ) tr CCT(9) A txc Hg) B (5.8) 

where the summation is over inequivalent (possibly disconnected) planar maps g, and n is 
the number of edges. 

Thirdly, we replace these maps g with their dual maps g. Note that this must be 
performed separately for each connected component. The correspondence is again one-to- 
one. An example is given on Fig. 2. The new maps have bicolored faces; equivalently, one 
can orient their edges so that white (resp. black) vertices correspond to clockwise (resp. 
counterclockwise) faces. 

The resulting maps g are very similar to the Feynman diagrams of sec. 5.1; however, 
they still lack the "dotted lines" which link together the various connected components. 



This is where Eq. ( |5.7|) comes in. Indeed to each term in the sum (|5.7| ) we associate one 
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particular set of dotted lines as follows: since IT' > n^wg), and the cycles of a(g)r(g) 
are associated to faces of g, one can think of IT' as ways of grouping together faces of g, 
that is vertices of g: these are precisely the dotted lines. Finally, the other conditions 
in the summation of Eq. ( |5.7| ) can be interpreted as follows: #(n' V V R T ( g )) = 1 

ensures that the diagram including dotted lines is connected; and #Il T ( 3 ) cr ( s ) — #n' = 
#(IL a (g) V II T ( S )) — 1 ensures that it is "minimally connected", i.e. that it has a tree 
structure. 

Finally, we compare the weights: the Wc p . are associated to groups of vertices of g 
linked together by dotted lines, just as in the previous section; as to the tic n(g) A (resp. 
trc r(s) -B), they associate to each white (resp. black) vertex of g, that is each clockwise 
(resp. counterclockise) face of g, a weight of 9 P (resp. 9 P ) where p is the number of edges 
surrounding it, which is again what is required. 

Remark: the special case Ha-(g) V Tl T (g) = HoV^Wg), for which the associated weight 
is a single Wc rj(g)T{g) , occurs when g is a disjoint union of bicolored trees, or equivalently 
when all vertices of g are linked together by dotted lines. 

5.3. Case of rectangular matrices 

We sketch here how the diagrammatic expansion of sec. 5.1 can be generalized to the case 
of the integral (|3.5| ) over rectangular matrices. Introduce complex rectangular N\ x N2 
matrices X and Y; then, applying the same trick as above, we find 

/ DU j DUexp — Tr{AUBV ] + VB ] U ] A ] ) 

JU(N 2 ) JU(N 1 ) s 

= const J DXDXtDyDyt e -^ Tr (^ t +^^ t )+^X^ t ^ t ^^)+^^( Bt ^ t ^ B ^) (5.9) 

where = min(A^i, N2), N' = max(iVi, N2). The diagrammatic expansion is identical to 
that of sec. 5.1, except that now there are two types of vertices, weighted by extra factors 
Ni/N' and N2/N 1 respectively. Due to the presence of X and Y', resp. Y and X\ in these 
vertices, we see that they must alternate and therefore these diagrams have both faces and 
vertices bicolored. Up to a factor of 2 corresponding to the two possible colorings of the 
vertices, this is the same as requiring diagrams to have bicolored even-sized faces. Each 
face of size 2p receives a weight Tr(AA^) p or Tr(BB^) p depending on its color (orientation). 

Note that if N± = N2 = N, the diagrams are weighted identically as in sec. 5.1; this 
means that the integral ( |5.2| ) for square matrices can be considered in the large iV limit 
as the particular case of the integral (|1.1|) for which all odd moments vanish (cf a similar 
observation in sec. 4.4). 
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6. Summary of results, Tables 

One may compute the first terms of the expansion F(A,B;s) = YunLi j^F n (9,9). It is 
sufficient to tabulate them for 9\ = 9\ = since if we write A = A' + 9± I, B = B' + 6\I, 
with A' and B' traceless, F(A, B; s) = F(A', B'\ s) + d\9i/s . Or alternatively 

F n (9i, 9i, 92, 92, ■ ■ •) = 9i9i8 n \ + F n (0, 0, 9' 2 , 9' 2 , ■ ■ •) 
with 9' p = J2 P q=Q Q)9 q (-9 1 ) p - q and likewise for 9 P . Up to order 8, we find: 

Fx=0 F2 = \ 02 02 i*3=3 03 03 

F4,=\[e 2 2 (30f-2e 4 )+6»4 (-zel+h)] 

ff>=i [e 2 9 3 (2oe 2 9 3 -58 5 )+e 5 (-5 2 3 +05)] 

f 6 =^ [e 2 3 (27 e|-i6 0J-3O §2 84+7 8 6 )+e 3 2 (-16 e|+6 S 2 +i5 e 2 4 -3 6 6 ) 

+3 2 4 (-10 0~f+5 0^ + 10 2 04-2 6 )+0 6 (7 0| -3 0^-6 2 04+0e)] 
F 7 =0 2 2 03 (66 02 03-21 2 05-20 03 4 +40 7 ) 
-03 04 (20 02 03-5 03 04-6 02 05+0?) 
-02 05 (21 B~l 3 -6 3 04-6 02 05+0V) 
+ i 07 (28 Ol 3 -7 3 04-7 02 05+0 7 ) 

F 8 = -| [3 2 4 (117 0^-192 02 0^-180 9\ 4 +25 0^+56 3 5 +56 2 6 -lO0 8 ) 
-4 2 3 2 (l44 0*-176 2 0^-200 0f. 4 +25 0^+48 3 5 +56 2 6 -90 8 ) 
-4 2 2 4 (135 0^-195 0^ 04+25 0^+54 03 5 +0 2 (-200 0^+56 6 ) -9 0g ) 

+ 4 2 (75 02-100 02 2 -lOO0 2 04 + 1004+24 03 05+28 02 6 -4 8 ) 

+0 3 5 (168 0^-192 02 0^-2160| 04+24 0^+48 03 5 +56 2 6 -80 8 ) 
+0 2 6 (168 0^-224 02 0^-224 0^ 4 +28 0^+56 3 5 +56 2 6 -80 8 ) 
+0 8 (-30 0^+36 02 03+36 0! 04-4 0^-8 3 5 -8 2 6 +0s)] 



One may slightly simplify this expansion by making the following asymmetric change 
of variables: we introduce instead of the moments 9 q , the free cumulants (p q which are 
defined by Hl^ig 



v (g + Ei^-2)! n (-fl 



ol^ , . . . ,OLq >0 



a,; 



(6.1) 



Indeed, if one now expands F as 



d 2 F 



2! ^ 89 a d9 r 

q,r>l q r 



0=0 



9q9 r + 
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then the first 3 derivatives were computed exactly in [BBJ using the formalism of dispersion- 
less Toda hierarchy, and turn out to be simply expressible in terms of the <j) q . In particular, 
§§-\ d _ = ~0<?- The same expansion to order 8 now takes the form: 

F 5 = jr9 5 5 -9 3 6 2 (4> 5 +4> 3 4> 2 ) 

^6 = ^606-0402 (06+04 02 + ^03) (06+04 02+03 + ^10 + ^2 (706 + 120402+503+202 1 ) 

F 7 =ie 7 <^ 7 -6l56l2(^7 + 05^2+^403)-e4e3(07 + 05^2+2^4^3+^302)+ e 3^ (40 7 +70 5 02 +804 03 +50302 ) 
-^8= 08— ^6^2 (08 + 06 02 +05 03 + 5 01)— #5^3 (08+0602+20503+01+04 02+03 02) 

-^l(0 g +0 6 02+20503+|01+0402+20302+i02) 

+ \B A d\ (908 + 160602 + 180503 + 1101 + 110402 + 140302+01) 

+6 2 3 d 2 (908 + 160602+240503 + 1101 + 16040^+200202+204) 

-|^(lO08+240 6 02+240503 + 150l+2404 0^+240202+30^) 

One notes the intriguing feature that in either basis, for all known terms (up to 
p = 10), pF p has only integer coefficients. In other words, the quantity sdF/ds has only 
integer coefficients. The rationale behind this fact and its reinterpretation in terms of the 
counting of "rooted" objects has remained elusive to us. The integral ( |1.1| ) still has a few 
mysteries. . . 

In this article, we have presented the standard lore on the integral ( |1 . 1| ) as well as 
recent developments. The integral itself is well understood by a variety of methods (sect. 2) 
and admits interesting generalizations (sect. 3). However the explicit expression, in terms 
of symmetric functions of the eigenvalues, of the integral itself (not to mention the more 
complicated problem of the associated correlation functions) is more subtle; sect. 4 was 
devoted to the discussion of this problem in the large iV limit. These considerations, based 
on integrable hierarchies, are complemented by other recent work based on combinatorial 
arguments, as discussed in sect. 5. It is our feeling that the exact interrelations between 
these various approaches require further study. 
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